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Abstract 

Two-photon absorption is theoretically analyzed within the semiclas- 
sical formalism of radiation-matter interaction. We consider an ensemble 
of inhomogeneously broadened three-level atoms subjected to the action 
of two counterpropagating fields of the same frequency. By concentrating 
in the limit of large detuning in one-photon transitions, we solve pertur- 
batively the Bloch equations in a non-usual way. In this way we derive 
an analytical expression for the width of the two-photon resonance that 
makes evident sub-Doppler two-photon spectroscopy. We also derive an 
analytical expression for the Stark shift of the two-photon resonance. 

PACS; 42.50.-p (Quantum Optics), 42.62.Fi (Laser Spectroscopy) 

Abstract 

Se analiza teoricamente la absorcion de dos fotones dentro del formal- 
ismo semiclasico de la interaccion entre la radiacion y la materia. Con- 
sideramos un conjunto, con ensanchamiento inhomogeneo, de atomos de 
tres niveles sometido a la accion de dos campos contrapropagantes de 
igual frequencia. Resolvemos perturbativamente las ecuaciones de Bloch 
del sistema de una forma no usual concentrandonos en el limite de alta 
desintonia de las transiciones a un foton. De esta forma obtenemos una 
expresion analitica para la anchura de la resonancia a dos fotones en la 
que se pone de manifiesto la posibilidad de espectroscopia sub-Doppler a 
dos fotones. Tambien obtenemos una expresion analitica para el desplaza- 
miento Stark de la resonancia a dos fotones. 

1 Introduction 

Two-photon absorption (TPA) is one of the most basic radiation-matter inter- 
action mechanisms. It consists in the excitation of an atom or molecule from 
a lower quantum state |1) to an excited state |2) of the same parity as |1) 
in a single step. In this case the initial and final states cannot be connected 



through an electric-dipole transition. Thus parity conservation imphes that two 
hght quanta must be absorbed simultaneously. The theory of TPA was first 
developed by Maria Goppert -Mayer in 1931 in her Ph.D. Thesis PP. 

As a multiphoton process, TPA is closely related to Raman scattering. In the 
latter process, one photon is absorbed while the other is simultaneously emit- 
ted, the energy difference being retained by the molecule. While spontaneous 
Raman scattering was observed as early as 1928 0, TPA was not observed until 
1961 13 after the advent of the laser (in fact TPA is one of the first nonlinear 
optical phenomena demonstrated with the aid of laser radiation). The reason 
for that delay in the observation of the two multiphoton processes lies in the 
fact that while in spontaneous Raman scattering the scattered light intensity 
is proportional to the intensity of the incoming radiation, in TPA the power 
absorbed is proportional to the square of the intensity of the incoming field and 
thus higher excitation energy is required for TPA. 

TPA is a very important tool in laser spectroscopy as it makes possible 
the transition between two states that cannot be connected by electric-dipole 
interaction. Of course these transitions can also be investigated by making use 
of resonant one-photon processes through an intermediate level, but in this case 
the measured linewidth of the process is increased by the linewidths of the two 
successive one-photon absorptions. TPA also allows the coherent excitation of 
molecules to states whose energies fall in the far ultraviolet, by making use of 
visible radiation, for which coherent sources are easily available. 

One of the most outstanding features of TPA is that it allows sub-Doppler 
precision measurements ^. This last fact was first analyzed by Vasilenko et al. 
0] in 1970 and observed in 1974 |S1E|- Doppler broadening comes from the 
fact that atoms moving with different velocities "see" the field with different 
frequencies because of the Doppler effect. This is a source of inhomogeneity that 
increases the measured absorption linewidth. In one-photon transitions this 
limitation cannot be easily overcome unless subtle phenomena such as the Lamb- 
dip produced by spectral hole burning are exploited. In TPA, however, there is a 
simple way of (almost) getting rid of Doppler broadening. This occurs when the 
two photons inducing the transition come from two counterpropagating beams 
of equal frequency. In this case all atoms are in resonance with the two-photon 
process since the Doppler frequency shifts of the two photons "seen" by the 
atom are opposite among them, independent of the atom's velocity. Hence 
the sum of the energies of the two photons, as "seen" by any atom, is twice 
the energy of a single photon in the laboratory frame, and the inhomogeneity 
almost disappears. 

In quantum optics textbooks, TPA is often introduced after field quanti- 
zation [7]- Nevertheless TPA does not need the existence of photons to be 
understood and some textbooks analyze the phenomenon from a semiclassical 
point of view [5j that is, by treating matter quantum-mechanically and radi- 

^ Raman scattering also allows the investigation of transitions in which the initial and final 
states are of the same parity. With respect to Doppler compensation, in Raman scattering it 
is only partial and the degree of compensation depends on the energy difference between the 
initial and final molecular states. 
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ation classically (in this scmiclassical approach one must understand that the 
word photon refers to the amount of energy absorbed by the atom, not to any 
quantum already existent in the free electromagnetic field). There are several 
ways of stiidying TPA in this scmiclassical approach: derivation of nonlinear 
susceptibilities, application of standard perturbation theory, even derivation of 
exact analytical results. Nevertheless to our knowledge no standard textbook 
derives the main characteristics of TPA (such as Dopplcr compensation and 
the Stark shift of the resonance) within the semiclassical frame. In this article 
we give a compact and clear presentation of TPA from a semiclassical point of 
view, by solving perturbatively the equations of motion for the density matrix 
elements. 

2 Semiclassical density matrix equations. 

Let us consider a classical monochromatic electromagnetic field of the form 

E {z, t) =e [El cos {(jt + kz) - E2 cos {ujt - kz)] (1) 

where e is the unit polarization vector (linear polarization is considered) and 
Ex and E2 are the constant real amplitudes of two counterpropagating plane 
waves of angular frequency u> and wavenumber fc, which travel along the z 
axis. Note that this form of writing the total field is completely general for 
the superposition of two counterpropagating monochromatic linearly polarized 
waves of equal polarization, since any dephasing between them can be removed 
by suitable choice of time and space origins. This field represents a standing 
wave when Ei = E2 and a traveling wave if either Ei or E2 is taken to be zero. 

This classical field interacts with a medium composed of three-level atoms 
(Fig. 1): levels |1) and |2) of the same parity, and contrary to that of the 
intermediate level |0). This is the simplest level scheme that allows the descrip- 
tion of TPA in terms of the usual electric-dipole Hamiltonian. In this way, the 

transition |1) < > |2) is produced via the virtual transitions |1) < > |0) and 

|0) < — > 1 2) (state |0) is kept far from resonance). The existence of an interme- 
diate level enhances the excitation probability between states |1) and |2) as will 
be shown. 

The unperturbed hamiltonian Hq of the three-level atoms is given by (see 
level diagram in Fig.l) 

Ho = n{u2o\2){2\-coi\l){l\), (2) 

and the origin of energies has been taken at the intermediate state |0). Since 
levels |2) and |1) have the same parity, and contrary to that of state |0), the 
allowed electric-dipole transitions are |1) ■* — > |0) and |0) < — > |2). Thus the 
dipole-moment operator is written as 

/X = /X20 |2) (0| + H02 |0) (2| + /xio |1) (0| + Hoi |0) (1| , (3) 
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where fi^^ = that can be taken to be real without loss of generality 

through proper choice of the basis states phases (/x^ = /i^j). The interac- 
tion hamiltonian of an atom located at z reads Hi {z,t) — — /i • E(z,i) and 
the total hamiltonian that governs the coherent evolution of the atoms is then 
Hs = Hq + Hi , where the subscript S is used to denote the Schrodinger picture 
implicitly adopted. Before solving the Schrodinger equation, it is convenient 
to remove fast oscillations at optical frequencies appearing in the hamiltonian. 
This is accomplished by transforming from the Schrodinger picture to the field- 
interaction picture. The appropriate unitary operator for making such transfor- 
mation is 

U (t) = e"^' |2) (2| + |0) (0| + e-'"* |1) (1| . (4) 

Note that this operator is similar to that defining the Dirac picture but, instead 
of removing the fast free atomic evolution (which would be accomplished with 
the operator Ud (t) = e"^^°* |2) (2| + |0) (0| + e-*"°i* |1) (1|), we remove the 
fast dynamics originating from the optical frequency of the field. In the new 
picture, the state vector \^) of the system (^\ip) = U {tps)^ obeys the following 
Schrodinger equation 

where the hamiltonian H in the new picture^ is calculated TU' through 

H = UHsU-^ +in—U-\ (5) 

After performing the rotating wave approximation [3 |3 (that consists in re- 
moving fast oscillating terms) the hamiltonian reads 

H (z) = h {-52 |2) (2| +5i\l) (1| - ^iE\2) (0| - E |0) (1| + h.c.) , (6) 

where 

E = cl^ie"'^ - <l>2e-''^ (7) 

Sl(2) = ^^-^^01(20)- (9) 

and h.c. stands for hcrmitian-conjugate. Note that the new picture, in combi- 
nation with the rotating wave approximation, yields a hamiltonian independent 
of time. E {^E) is half the complex Rabi frequency of the field associated with 
the lower (upper) transition of an atom located at position z. 

Now we determine the evolution equation of the density matrix (more prop- 
erly: the population matrix, see below). We choose to work with the density 
matrix instead of the state vector since in this way relaxation and pumping 

■^Under a unitary transformation, any operator Os (in the Schrodinger picture) transforms 
according to the rule O = UOsU~^ . Notice that H is not transformed in the same way. 
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processes can be (phenomenologically) incorporated into the model in a sim- 
ple way. As we are considering not a single molecule but a large number of 
molecules which are moving at different velocities, an ensemble average must be 
performed. The ensemble averaged density matrix is usually called population 
matrix j^]. This ensemble must be defined for each velocity and, since the 
interaction depends on space, the population matrix must also be defined as a 
function of the position z: 

p{v;z,t)^U{vy^Y.Pa{^-:Z,t)- (10) 

a 

Here p is the population matrix, is the density matrix for an atom labeled 
by a, and a runs along all molecules with velocity v that, at time t, are within 
z and z + dz . J\f (v) is the number of such molecules, which is assumed to be 
independent of z and t (homogeneity and stationarity of the velocity distribution 
is assumed) . The equation of evolution of the population matrix is formally like 
the Schrodinger-von Neumann equation governing the evolution of the density 
matrix of a single atom, plus an additional term 

{dt+vd,)p,^^{th)-^\HM +(tp) , (11) 

= 0,1,2). Fp describes irreversible processes (relaxations and pumping) 
and r is a generalized Liouvillian. In this article we shall consider the simple 
expression 

(fp)^^ =-7P,^.+7<5,,i5,,i, (12) 

with 5 the Kronecker delta. The first contribution describes relaxations in a sit- 
uation in which all density matrix elements decay with the same constant 7 (ab- 
sence of dephasing collisions jllp. The second contribution (pump) guarantees 
that the ground state |1) is asymptotically filled in the absence of interaction. 
With this choice for Fp, Tr (p) = 1 always. We adopt this simple limit because 
the expressions are much clearer and the details of the relaxation processes do 
not modify the essential physics of TPA. 

By substituting Eqs.© and (|12|l into Eq.l^TJ, the final equations of evolution 
of the population matrix elements run 

{dt+vd,_)p22 = -iP22 + H^{Epm- E*p^„) , (13) 

((9f -|-w9^)poo = --/Poo + iiEpio-E*Poi)-ip{EpQ2-E*p2Q), (14) 

{dt+vd,)p,, = ^{\- p,,)+i{E*p,,-Ep,,), (15) 

{dt+vd,^)p2i = -(7-i(5)p2i +i£'(ppoi -P20), (16) 

{dt+vdz)p2o = -[l-i^-i^^ P2n+ipE{pQ^- p22)-iE*p.2^, {11) 

{dt+vd,)pai = - (1 -i^AAl Poi+iEipii- Poo) + itJ-E*P2i- (IS) 
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where 



A = 61- S2 = UJ20 - t^oi, 



(19) 
(20) 



have the meaning of two-photon detuning and intermediate level detuning, re- 
spectively (see Fig.l). The above equations should be complemented with the 
evolution equation of the electromagnetic field. Nevertheless we shall consider 
-E as a parameter. This corresponds to a physical situation in which the gas of 
molecules is confined within a small region of the space which is large compared 
with the radiation wavelength but small enough for neglecting field depletion 
(thin film approximation). 

Note that A is a structural parameter of the atoms, and we shall consider 
only the case in which A is a very large quantity as compared with the rest of 
frequencies (7, S, E) appearing in the problem. This limit guarantees that one- 
photon processes (i.e. the electric-dipole transitions |1) < — > |0) and |0) < — > |2)) 
are severely punished since the one-photon detunings 5i (« A) and 82 (~ —A) 
are much larger than the widths of the one-photon resonances. For example, 
consider the states |2) = 851/2, |0) = 7P and |1) = 65*1/2 of Cs. In this case |T2] 
woi = 4.098 • lO^^s-i and LO20 = 0.489 • W^s~^ and thus A = -3.609 • W^s'^. 
Cs is a gas and the one-photon transitions width can be estimated to be given by 
their Doppler width which, at room temperature are (see Section 4) 2.6 • 10®s~^ 
and 2.21 • lO^s"^ for the upper and lower transitions, respectively: in this case 
there are nine orders of magnitude between A and the width of the one-photon 
resonances. 

Eas. (|13|l - H18|) do not admit a simple analytical solution in the general case 
but can be solved perturbatively in the case of very large A. This is done in 
the next section. 

3 Perturbative solution of the steady state 

In this section we solve perturbatively the equations of evolution of the density 
matrix in steady state {dt — > 0). Note that this is the state asymptotically 
reached by the system due to the presence of relaxations. Here we present 
the main results and leave the details to Appendix A. As commented, we shall 
consider the limit A 3> 7, (5, E. We also consider that the inhomogeneous width 
7„ (see Section 4) is small as compared with A, that is, we assume that kv 
with k the field wavenumber. This can be made formally explicit by writing 
A = e~^Ai with Ai a quantity of the same order of magnitude as the rest of 
the frequencies present in the problem and < e <C 1 (smallness parameter). 
We also make a series expansion of the density matrix elements of the form 



00 




(21) 
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Substituting this expansion in the population matrix equations and identiiying 
terms of equal powers in e, one gets 

{vd.+l)p^^-i{Ep'(^-E*p'^^) (23) 

-7 + {vd. + 7) P^i^ + ^ {ep^^ - E*p^^) , (24) 
+j-i5) pif - iE (pp^^ - p^^) , (25) 
{vd, + 7 - i 1 5) - ipE - p^^) + iE*p^^'> , (26) 
{vd^ + 7 - i 1 5) - iE - - ipE*p^^^ , (27) 

1 to oo. Note that these equations refer to an ensemble of 
atoms moving with velocity v located at z. These equations can be solved at each 
order n of s. Wc can integrate the first four equations to obtain the populations 
pI"'' and the two photon coherence pj"'' if we know the value of the one photon 
coherences at this order {pqi and P20 ) • These quantities are obtained from the 
two last equations. Note that the form of these last two equations (which relate 
two consecutive orders) allow the values of the one-photon coherences at a given 
order n + 1 to be algebraically determined in terms of the previous order n. In 
particular, for n = — 1 we obtain Pq°^ = pjo' = 0, since PiJ^^ = 0. These values 
allow to solve, from the first four equations, for the rest of matrix elements at 
order n = 0. Next, p^^ and P20 determined from the last two equations and 
so on. There is just a single point that deserves some explanation and concerns 
the integration in z of the first four equations. Notice that although wc do not 
know any boundary conditions (in terms of z) for the variables, we can make use 
of the knowledge that, when the field is off {E = 0), all variables must vanish at 
any order but p^^ , which must be equal to unity since Tr {p) = 1. In Appendix 
A the equations are solved systematically. In the following we make use of the 
result of the integration. 

4 Velocity and space averages 

We must concentrate on the calculation of a quantity directly related with mea- 
surement. We shall consider the fluorescence intensity from the system, which 
is directly proportional to the amount of population excited to the upper level. 
The fluorescence signal collected by a detector will come from all atoms (all 
velocities) existing within a finite region (of length L) of the system. Thus it 
is necessary to perform both spatial and velocity averages. The spatial average 



= 
= 

= 
= 

where n runs from — 
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reads 



(28) 



where we shall take, as already commented, L ^ X (with A the light wavelength) 
since typically the detector will collect the fluorescence from a "macroscopic" re- 
gion of the system. It is evident that only the spatial dc component of P22'' ("^^j z) 
will contribute to the spatial average since L 3> A, as stated. Consequently 
it will suffice to calculate only those terms. 

With respect to the velocity average a few words are in order. In a gas, 
inhomogeneous broadening is due to the Doppler effect which is different for each 
atomic velocity. The atomic velocities of a gas obey the Maxwell-Boltzmann 
distribution 



g{v) 



1 



: exp 



2k /In 2 

Iv 



exp 



2kv 



(29) 



with u the most probable velocity given by it = y^fceT/m {ks is Boltzmann's 
constant, T is the absolute temperature, and m is the molecular mass). 7^ = 
2V^n 2ku is the inhomogeneous HWHM (half-width at half maximum) in terms 
of the frequency f2 = 2kv (the factor 2 is added for later convenience, since in 
TPA it is not the radiation frequency -or its wavenumber- that is the important 
parameter but twice its value). The problem with the Gaussian distribution 
is that some integrals appearing in the final expressions cannot be evaluated 
analytically. In order to obtain analytic expressions as simple as possible, we 
shall consider a Lorentzian distribution for the atomic velocities 



1 



Civ) = - 



2k 



Iv 



{2kvY 



(30) 



where 7^ = 2ku is the inhomogeneous HWHM (half-width at half maximum) 
in terms of the frequency SI = 2kv. The results obtained with this distribution 
will differ quantitatively but not qualitatively from the Gaussian distribution, 
as will be shown. 

The averaged population of the excited level is then calculated through 



P22 ) = 



dv 



(") r \ 



C{v) 



(31) 



Clearly the averaging order is unimportant. We could first perform the velocity 
average and then the spatial average, obtaining the same result. From the 
computational viewpoint however it is more convenient to perform first the 
spatial average since in this way the ac-components (in terms of z) of are 
removed from the calculations from the beginning. 

From Eas. (|84|) and (|86|l of Appendix B, the fully averaged population of the 
upper level reads, up to order e'^. 



3 / J3)\ _ 
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where 



and 



No 



x2 \ 3 



1 / 1 



{.ay 



(l + 7j' + f 



In writing Eas. (|32|) and (|33|l we have introduced the notation 
and the normahzed frequencies 

1 1 



(32) 



(33) 



(34) 



(35) 



Note that e" combines with Aj"" in both orders n = 2 and n = 3 to yield 
A~", leading to a final expression independent of e. Next we analyze these 
expressions. 



5 Analysis of the results 

5.1 Strength and width of the resonance 

In order to analyze the strength and width of the resonance it is suffice to 
consider the dominant term 7V2. General results are: (i) TPA is proportional 

to the squared field intensity (f/*^)^, (ii) The existence of an intermediate level 
with a finite detuning A enhances the probability of the process (the smaller 
A the larger amount of excited population) , and (iii) The maximum transfer of 
population is produced at (5 = (this result will be corrected at the next order; 
see next subsection). Let us consider some special cases. 
In the case of homogeneous broadening (7^ = 0), 



+ 4^2 + 1 



(36) 
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Note that N2°™ is proportional to {A^ + AA^ + 1), which in its turn is propor- 
tional to the mean value of the squared field intensity -a signature of two-photon 
absorption. This factor is six times larger for standing waves than for travel- 
ing waves. Note that this (important) numerical factor is the single difference 
between standing and traveling wave configurations in this homogeneous broad- 
ening limit. We conclude that, from an experimental point of view, it is most 
convenient to illuminate the cell with a traveling wave and make it reflect on a 
mirror located after the cell in order to produce a standing wave. This repre- 
sents no extra energetic cost and the fluorescence signal collected in this way is 
6 times larger than without the mirror. 

With a non-null inhomogeneous broadening two limits of interest are: a) 
excitation with a traveling wave {A = 0) 



(1 + 7J 



and b) excitation with a standing wave (A — 1) 



7A 



2(1 + 7.) 
+ + f 



(37) 



(38) 



Note that the effect of the inhomogeneous broadening is dramatically different 
for traveling wave or for standing wave cases: if 7^, 3> 1 {i.e. 3> 7, inhomo- 
geneous limit) TVj^ 0, whereas in the same limit ^^2^^^ ^N^°'^'^^ , where 



N. 



horn, SW 



is obtained from N^""^ 



by putting A — 1. 



In order to make clearer comparisons among different cases we next analyze 
next the maximum of iV2 (that occurs at J = as stated) and its width in terms 
of 5. 

The maximum of N2 reads 

2 



^max 



7A 



(A4 + 4^2 + 1) + 47,^2 



(39) 



(1 + 7.) 

In Fig. 2 N^'^^ (normalized to its maximum value, for A = 1 and 7^ = 0) 
is plotted as a function of the inhomogeneous-to-homogeneous widths ratio 7^ 
for different values of A. Clearly, for A = 1 (standing wave) TPA is almost 
insensitive to the amount of inhomogeneous broadening, whereas ior A — 
(travelling wave) the decrease in TPA is dramatic for ratios as moderate as 
7^, = 1 or larger. 

As a function of the normalized detuning 6, N2 has a bell shape whose 
FWHM r is easily calculated from Ea. (p!^ . and reads 



p2 



^w+{w-lf p + {w-l)f 



(l+7.„)': 

\ \ + A^-A{l + ^,] 

2 1 + ^4 + 4 (1 + ^ J ^2 ■ 



/ = 



A^ 



(40) 
(41) 
(42) 
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For a homogeneously broadened line (7^ = : w = 1) the width reads Fhom = 2 
{i.e. in terms of the detuning S the width reads 27). In the special case of 
a traveling wave (A = : / 1/2) the width reads V'^^ = 2 (1 + 7^), i.e., the 
width is the sum of the homogeneous and inhomogeneous widths. For a standing 
wave {A = 1) no simple expression is obtained; nevertheless in the special case of 
large inhomogeneous broadening (7^ ^1:/^— 1/2, 00) the width reads 
F'^'^ (7„ — > 00) = 2 (1 + l/27„) which tends to the homogeneous width 2 for 
sufficiently large inhomogeneous broadening. This last result is a fundamental 
property of TPA: sub-Doppler spectroscopy can be performed in TPA experi- 
ments by using a standing wave 00 E El IS. 9.- 

In Fig. 3 we plot F/homj as given by Ea. (|40ll . as a function of the normalized 
inhomogeneous width 7^ for ^ = 1 and A = 0.5. Clearly, for any value of A 
different from zero, F/hom first grows until it reaches a maximum and finally 
decreases tending to unity for large enough 7^. Of course the optimum situation 
corresponds to ^ = 1'^. Thus for large enough 7,y the inhomogeneous broadening 
does not contribute at all to the width of the resonance. 

In Fig. 4 we show the same representation for A = 1 (full line) together with 
the numerical integration assuming a Gaussian velocity distribution. It can be 
seen that the dependence is qualitatively the same and that only relatively small 
numerical deviations are appreciated between both cases. This confirms that 
the exact form of the velocity distribution is not very important, whenever it is 
bell shaped. 

5.2 Shift of the resonance 

As we have seen, at order the maximum of the resonance is located at S = 
0. Nevertheless, two-photon processes induce a shift of the resonance, the so 
called Stark shift. This shift is only captured at third order of the perturbative 
expansion. Making use of Eqs.(j22|) and H33|l . we compute d {N2 + N3) /d6 = 
and obtain 



.l^J (1 + A4) + 4A2 (1 + 7j3 

(43) 

which is the Stark shift. Note that this shift is proportional to (0^/7A), and is 
thus of order e. Note also that whenever /i = 1 (i.e. when both one-photon tran- 
sitions have equal electric dipole matrix elements, see Eq.®) the shift vanishes. 
We see that the sign of the shift depends both on the sign of the intermediate 
level detuning A and on the asymmetry between both one-photon transitions 
through the quantity (/J,^ — l) . 

^ A 7^ 1 can bo understood as the sum of a traveling wave and a standing wave. Thus 
the result in that case is the sum of the two contributions. As the TW contribution is less 
important the larger is 7,, and the contribution of the SW is basically independent of 7^ this 
explains the above result. The main difference between A = 1 and A I lies in the strength 
of the resonance as shown in Fig. 2. 
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Particular cases of interest are: a) excitation with a traveling wave {A = 0) 

sZr, = 2 (m^ - 1) (;^) , (44) 

which is independent of the inhomogeneous broadening, and b) excitation with 
a standing wave {A = 1) 

which tends to Sstarkl'^ for 7„ ^ oo. In Fig.5 Sstark/htark is represented as a 
function of the inhomogeneous width for both Lorentzian broadening (Ea. H43|) l 
and Gaussian broadening. Again it can be appreciated that the results are very 
similar for both types of inhomogeneous broadening. 



5 + 37^ + T-g 
l + 2(l + 7j' 



(45) 



6 Conclusion 



In this article we have analytically studied two-photon absorption (TPA) in 
an inhomogeneously broadened medium pumped by two counterpropagating 
light beams of equal frequency. By making use of perturbative techniques, we 
have derived explicit analytical expressions for the strength and width of the 
resonance as well as for the Stark shift in the case of Lorentzian broadening. 
Comparison with Gaussian broadening (numerically computed) has shown that 
the qualitative features of TPA are quite independent of the specific type of 
inhomogeneous broadening. 



7 Appendix A 

At order one trivially gets 

pf} = pf^ = 0. (46) 
At order e'^ the equations are 

vdj^,^ = ~jpg\ (47) 

vdJ^J - -ipit (48) 

vdj^,^ = 1-7p1i\ (49) 

= (50) 

vdj^^ = -^Api]^+^E^^[plS-pg^)~^E*p^^,\ (51) 

vd.pl^^ = +'^^pl^^'+^E{p^|^-p^^^)-^E*pp^^,\ (52) 
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whose solution is 



Pii V", z) 

(0) ( \ 



Poo ^) = P21 = 0> 



and 



Poi ^) 

P20 ^) = 0- 
At order the equations are 



vdzPoo 

vdzP2o 
vdzPoi 



= -7P^o' + * ( Ep'^^ - E*p\;^ ] - ifi ( E^ip);^ - E*p)^^ 



.(1) 



(1) 



.(1) 



,(1) 



.(1) 



.(1) 



(53) 
(54) 

(55) 
(56) 

(57) 
(58) 
(59) 
(60) 



= - ( 7 - ) P^o^ - ^ AiP^o^ + ^Ep ( p'^> - p'^> ) - iE*(^ ,(61) 



.(2) 



.(1) 



.(1) 



P^V + -.^iP'ol + (piV - P^o^ ) - *i^>P^V,(62) 



.(2) 



.(1) Jl) 



and integration along z has to be carried out. By using Eas.H55|l. it is straight- 
forward to obtain that 



(1) ( \ 

P21 K"^^) 



0, i = 0,1,2 
lip 



Ai 



and 



(2) / \ 

P20 yv, z) 



(2) ^ \ 

Poi l^'' 2:) 



with 



4i/i '/'?<?!>2„3jfcz , '/'i , 20102 \„jfc^ 



6^ 



4»p' r 
A? ^ 



20102 



Do 

a2 



„~ikz I '^l'^2 -3»fcz] 

(7 + i?+)0i , 



(7 + ^-)'/'2 

2/x2 



2^2 

20102 
^0 



01 ^ 20102 
^0 



^4 
ikz 



nV2_^~3ikz 



D. 



D± = 7 - i ((5 T 2/ct;) , 
= ^ — iS. 



(63) 

(64) 

(65) 

^rkz 

X,66) 



(67) 
(68) 
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At order jt is not necessary to compute all the terms since we are only 



interested in p^^2 ^-nd P20'' (^'^^ latter is necessary for calculating p)^^ at order 



.(3) 



J3) 



The necessary equations are 

a (2) (2) , ■ I (2) (2; 

vd:,p\{ = --fpy + I p [Ep'J - E p^o' 



a (2) 

(2) 

a (2) 



(2) 

- 7Poo 



7- 



P20 



.Ai (3) 



(2) 
02 



(2) 
P20 



^P-E I Poo ^ P22 



(69) 
(70) 
(71) 



-zi?*p('i',(72) 



and the searched quantities are given by 
4/i2 



(2) ^ \ 
P22 i'"^ Z) 



2 j,2n 

1<P2 



7A? 

16^2 



'2 / 7 + 



A? 



2ikv 



+ C.C. + 



(73) 



DID+ DoD*_ 



i2kz 



+ C.C. + 



-terms with e^"'^'"'' 



(2) / \ ^ 

Poo i^^ 2:) = -2 



A? 



7 + ikv 



i2kz 



7 



(2) / \ 

P21 iv, z) 



+terms with e*"*"^ 

4^0102 

AfDo 



2ikv^'^^ 
(n^0,±2) 



c.c. 



7 + + (m' - 1) 2 



2m'/'i 
AlDo 



7 + ^^ 



7 + £>_ + 2 - 1) 



01 + </'2 






^0 


^+ 












01+-? 




V ^0 


















V ^0 







(74) 



-(75) 



A2kz 



i2kz 



+terms with e'"'"'' {n^0,±2) 
from the three first equations and, from the last equation, 

pg^ {v, z) = p^o^^e*'^-^ + pgf\-'^'- + terms with e™'=^ 
where 

\p'-i) V 



n 7^ ±1 



(76) 



A? 



Dl 



\D, 



24 



+ 


'{p'~l) (23 


Do\ 


4p2 


/ 4 (7 + 


ikv) 




Do \Do D+ 


D\) 


^0 


\Do D+i^ 


- 2ikv) 



^102 



^Notice that if one is interested only in the analysis of the strength and width of the 
resonance (and not of the Stark shift), it is enough to calculate the non-oscillating term in 
Eg. 1731 by direct substitution of I65i in Eg. 1691 . quite a simple task. The rest of the terms 

(3) 

are necessary for obtaning of which becomes a simple but tedious task. 
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and 



(3,- 
P20 



" A? ■ 



1) 



+ 



4- 



Do 



2 



2 



3 



1 



202 



1?^ 



(7 + ikv) \ 
Do {j-2ikv)) 



4}, 



4 

D~o 



(7 — ikv) 
£»_ (7 + 2i/cw) 



7 



£1, 







7 — 2ikv 

1 2 

D^^D'o 



131 



1 



(7 — ikv) 
Do (7 + 2ikv) 



4}- 



(77) 



(78) 



Finally, at order we are only interested in obtaining the value of the 
population of the excited level. Thus we only need 



(3) 



(3) 
^ 7P22 

(3) 



iKEpo2 



(79) 



The spatial dc-component of P22 finally reads 



(3) 

P22,dc 



A3 

5 — kv 



2"<^l 



V l^+l 

5 + kv 



(01 + ^2) 01 



(80) 



+ 



( 5 + 2kv\ 



2"<^2 



(01 + 02) 02] 



\Do\ 



*102 



\Do\ 



8 Appendix B 



At order the spatially-averaged population of the excited level, Ea. (|73|l . is 

4^2 



7A? 



D, 



+ c.c. 



(81) 



Now the averaging over velocities has to be carried out. As v only appears in 
p'22 through D± (v) the only integrals to be done are of the type 



inti 



1 



d{2kv) 



Iv 



7 



72 + (2/fcw) 72 + (,5 ± 2/cw)' 



(82) 
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whose result is 



int — 



1 + lv 



and thus the averaged upper level population results to be 



(2)\ S^i 



1 + lv 



■2j2 1 



(83) 



(84) 



At order the situation is similar. Now the integrals that appear when 
making the velocity averaging of Ea. H80(l are of the type inti and also of the 
type 



1 



int2 (n) = — d(2kv) 



Iv 



{2kv) 



72 + ((5 ± 2kvY 



(85) 



(n = 1,2) whose result is 
int2 (1) = 

mt2 (2) = 
The final result reads 



lv5 



Iv^ (7 + 7^,) (37 + Iv) + ^'^ 



272 



16m' (1 + A^) ip" 



1)5 



2 /l2 



V7A1 
1 



+^ 

7J(7 + 7J'+^1 l^ol 



(86) 



27^(1 + ^4) 



(7 + 7.)' + '5' 
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Figure Captions 

Fig.l. Energy level diagram of the three-level atoms considered in the model. 
See text. 

Fig. 2. Maximum value of the population excited to the upper atomic level 
as a function of the inhomogeneous to homogeneous width ratio 7^/7 for several 
values of A. TV™"*^ is normalized to its maximum value (that corresponds 

to a homogeneously broadened medium (7^ = 0) pumped by a standing wave 
(A=l)). 

Fig. 3. Width of the two photon resonance normalized to the homogeneous 
width as a function of 7^/7 for two values of A. (Notice that for a traveling 
wave, A = 0, the width grows linearly with the inhomogeneous width as 27^/7. 

Fig. 4. As in Fig. 3 for a standing wave {A = 1) for both Lorentzian (full line) 
and Gaussian (dashed line) velocity distributions. 

Fig. 5. Dependence of the Stark shift on 7^/7 for a standing wave {A = 1) 
for both Lorentzian (full line) and Gaussian (dashed line) velocity distributions. 
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"Two Photon Absorption'', Perez- Arjona et al. 
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